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Abstract. We derive a recursion relation in the framework of Lagrangian perturbation 
theory, appropriate for studying the inhomogeneities of the large scale structure of the uni- 
verse. We use the fact that the perturbative expansion of the matter density contrast is 
in one-to-one correspondence with standard perturbation theory (SPT) at any order. This 
correspondence has been recently shown to be valid up to fourth order for a non-relativistic, 
irrotational and dust-like component. Assuming it to be valid at arbitrary (higher) order, 
we express the Lagrangian displacement field in terms of the perturbative kernels of SPT, 
which are itself given by their own and well-known recursion relation. We argue that the La- 
grangian solution always contains more non-linear information in comparison with the SPT 
solution, (mainly) if the non-perturbative density contrast is restored after the displacement 
field is obtained. 
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1 Introduction 

Analytic techniques for studying the inhomogeneities of the large scale structure (LSS) usually 
rely on several assumptions and approximations: (a) the LSS is formed due to the evolution 
of gravitational instability only; (b) the equations of motion are solved for a pressureless 
component of cold dark matter particles in terms of a perturbative series on an exactly 
homogeneous and isotropic background; (c) the use of the Newtonian limit, i.e., we demand a 
non-relativistic fluid, restrict to subhorizon scales, and assume negligible curvature; (d) there 
are no vorticities on sufficient large scales, and also primordial vorticity is absent; (e) the 
use of the single-stream approximation, i.e., neglecting velocity dispersion and higher-order 
moments of the distribution function; (f) the smoothing volume over the spiky Klimontovich 
number density is set to zero, i.e., neglecting backreaction effects on the velocity dispersion 
and on the gravitational field strength. We shall consider these restrictions in the current 
paper. They are appropriate for studying the weakly non-linear regime of structure formation. 
For departures of this framework see for example [1-8]. 

Perhaps the most straightforward analytic technique is called Eulerian (or standard) 
perturbation theory (SPT), since the equations are evaluated as a function of Eulerian coor- 
dinates [9]. Here, the local density contrast 5(x,t) = [p(x, t) — ~p(t)]/~p(t) and the velocity field 
of the fluid particle are the perturbed quantities. Importantly, the series in SPT relies on 
the smallness of these fields, and therefore breaks down as soon as their local values deviate 
significantly from its mean values. 

A convenient way to circumvent this drawback is to use the Lagrangian perturbation 
theory (LPT) [10-18]. In LPT, there is only one perturbed quantity, namely the displacement 
field It parametrises the gravitationally induced deviation of the particle trajectory field 
from the homogeneous background expansion. Therefore, the LPT series does not rely on 
the smallness of the density and velocity fields, but on the smallness of the deviation of the 
trajectory field. Furthermore, the explicit extrapolation of the Lagrangian solution leads to 
improved predictions even in the highly non-linear regime, whereas the series in SPT fails by 
construction [19]. 
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Historically, a great advantage of SPT with respect to LPT was the discovery of a simple 
recursion relation [20], whereas it was widely believed that there is no recursion relation in 
LPT [9]. A common argument for this absence was the additional complicacy in LPT, that 
Lagrangian transverse fields are also needed to provide an irrotational motion in the Eulerian 
frame. In this paper we derive an easy expression to maintain the Eulerian irrotationality, 

(n) 

thus constraining the Lagrangian transverse fields at each order n. Furthermore, we 

in) 

derive a relation to constrain the Lagrangian longitudinal fields ^ L , such that finally one can 
construct the displacement field \l/ = ^ n \I>( n ) in terms of the aforementioned longitudinal 
and transverse fields: * (n) = *[ n) + *^ n) . 

The LPT recursion relation is based on the fact that the nth order density contrast 5^ 
in SPT and LPT are in one-to-one correspondence with each other (while restricting to the 
initial position limit) [21]. To obtain the displacement field \I>( n ) we shall Taylor expand the 
Lagrangian mass conservation 5 = 1/J(*&) — 1, where J is the Jacobian of the transformation 
from Eulerian to Lagrangian coordinates. Importantly, by performing the Taylor expansion 
one loses the power of the non-perturbative formula. However, this is just a calculational 
step to obtain the displacement field — due to the inherent non-linearity in 1/J(\1/), the final 
use of the Lagrangian result should be concentrated on the unexpanded <5-relation. This 
point was first noted by Zel'dovich [19] who obtained an approximate solution by explicit 
extrapolation far into the non-linear regime. 



2 Formalism 



In the Lagrangian framework the only dynamical variable is the displacement field It 
maps the fluid element from its initial Lagrangian coordinate q to the Eulerian coordinate x 
at cosmic time t: 



x(q,t)=q+¥(q,t). 



(2.1) 



We utilise the Jacobian J = det[dx/dq] to describe mass conservation for our non-relativistic 
fluid: 



d 3 x = J(q, t) d 3 q , p(x, t) d 3 x = p(t) d 3 q , 



(2.2) 



where p(x,t) is the Eulerian density field, and p(t) is the mean mass density. Defining the 
density contrast 5 as 



p(x,t) = p(t)[l + 5(x,t)] , 



(2.3) 



we can use the mass conservation (2.2) to relate the density contrast to the displacement 
field: 



5(x,t) 



1 



J(q,t) 



1. 



In Fourier space this is [18] 



S(k,t) = / d 3 qe ikq [ e ife -*(9.*) - 1 



(2.4) 



(2.5) 



where we have used eqs. (2.1-2.4) for the last equality. We shall use the above equation below 
to relate the series in LPT to its counterpart in SPT. 
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In LPT, \I/ is expanded by a perturbative series, composed of purely longitudinal per- 

~ (n) ~ (n) 

turbations, labeled i 5! L , and purely transverse perturbations, denoted ^ T . As mentioned 
before we require an irrotational motion in the Eulerian frame, but Lagrangian transverse 
fields are mandatory to maintain this constraint. It is important to note that there is no 
decoupling between the transverse and longitudinal components. Thus, transverse and lon- 
gitudinal fields depend on both transverse and longitudinal fields. In Fourier space, the 
perturbation ansatz is: 

oo oo 

*(fc,t) = ^* (n) (fc, i ) = ^[*i n) (fc,i) + *S" ) (fc,t)] • (2.6) 

n=l n=l 

The solution for the fastest growing mode is 

^\k,t)=-iD n (t)L {n \k), (2.7) 
*5? } (fc,i) =-\D n {t) T (n) (fc), (2.8) 

where D = a is the linear growth function in an Einstein-de Sitter (EdS) universe normalised 
to unity at the present time to, and 

i3^_ . . . j3„ 



L (n \k) = J d Pl {2n) sn PU (2vr) 3 4 3 V-n - k) gWfo, ■ ■ ■ ,p B ) ~6 (Pi) • • • frfo.) , (2.9) 

T (n \k) = f dipi {2 ^n Pn (2vr) 3 ^g ) (p 1 ... w - k) SP( PI , . . . ,p n ) ~5 (pi) ■ ■ ■ 6o( Pn ) , (2-10) 

where we have employed the shorthand notation P\... n =Pi+P2+' ' '+Pn> an d Sq = S^ (t=to) is 
the linear density contrast; the vectors and are the symmetrised longitudinal and 
transverse kernels respectively, and they reflect the mode-couplings induced by non-linear 
evolution. The case n = l denotes the Zel'dovich approximation [19], i.e., we restrict to the 
initial position limit [21], and in this case the perturbative kernels are simply S^\p)=p/p 2 
and Srp (p) = 0, with \p\ = p. The Zel'dovich solution for the Fourier transform of the 
displacement field is thus 

* ( l\k,t) =-iD(t)L (1) (k) = -iD(t)S%\k)8o(k) = -iD(t)^~5 Q (k), (2.11) 
^x\k,t) = 0, (2.12) 

where we have used eqs. (2.7)-(2.10). In general, we expect longitudinal vectors of the form 
S™ = Pi2—n/Pi2—n^L^\ where is a scalar function to be determined, and the (non-zero) 

transverse kernels are constrained by P\2... n • = 0. In the following, we describe how to 

~ (n) (n) 

obtain solutions for the longitudinal and transverse displacements, \l/ L and ^ T for n > 2. 
3 From SPT to LPT 

In SPT, the density contrast is a perturbed quantity, and the solution for an EdS universe 



is: 



~5(k,t) = J2D n (t)~5^(k), (3.1) 
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with 

t (n) (k) = J d3pi {2 ^3 d n Pn W^iPl-n ~ fc)4 g) (Pl, • • • ,Pn) &(Pl) ' ' ' *>(Pn) , (3-2) 

where the symmetrised kernels Fn^ are given by the well-known SPT recursion relation [9, 
20, 22]. But note the special case f[ s) = 1. Now, using the above results for the displacement 
field and the density contrast, we Taylor expand eq. (2.5) and collect all individual nth order 
terms in eq. (2.5). The nth order density contrast reads then 

S (n) (k) = J d3 ^ )3 f ~ (2vr) 3 ^ 3) ( Pl ... w - k) XjMGft, . . , >p j ~ So(pi) . . . ~ So{pn) , (3.3) 

where we have defined the symmetric scalars Xn , which can be found in [21] up to fourth 
order. Explicitly, for n = l,2 and 3 they are 

x[ s \ Pl ) = k-S { ll T ( Pl ), ( 3 - 4 ) 
xi s \ Pl ,p 2 ) = k ■ S ( V T ( PllP2 ) + \k ■ S { H T ( Pl ) k ■ S% T ( P2 ) , (3.5) 

X 3 ] (Pl,P2,P3) = k ■ S lIt(Pi>P2,P3) + \ k ■ S lIt(Pi) k ' ^StCPs) k ' ^©tCPs) 

+ \ { k ■ S lIt(Pi) k ■ S l®t(P2,P3) + tw ° perms.} , (3.6) 

where we have introduced the short-hand notation = S^/ + S^. The Dirac-delta in 

eq. (3.3) fixes k = P\... n and as an important note we recognise the equivalence of the density 
contrast in both SPT and LPT: 

X«( Pl ,...,p n ) =F r [ s \ Pl ,..., Pn ). (3.7) 

«=Pl -n 

Above relation has been proven to be valid at least up to the fourth order [21]. Physically, 
this result is not surprising, because the transformation to Lagrangian coordinates should 
not alter the density contrast (neither the velocity of the fluid particle changes) . It is thus a 
reliable assumption to demand eq. (3.7) to arbitrary order, which we do so in the following. 

In our setting, the requirement of an irrotational fluid motion in Eulerian space implies 
vanishing (Lagrangian) transverse contributions up to the second order, i.e., S^L T = 

(2) (2) 

and S L( ^ T = S L . We will discuss this issue in the following section. This means that we can 

immediately obtain the longitudinal solution k ■ by equating eq. (3.5) and F 2 with the 
constraint k= Pl2 and the use of the Zel'dovich approximation, i.e., eqs. (2.11) and (2.12): 

Pn • sf{ Pl , P2 ) = if) (Pl ,p 2 ) - \ (l + ^) (l + =: B?\ Pl , P2 ) , (3.8) 

and as the s originate from longitudinal perturbations (i.e., the field perturbations can 
be written in terms of a potential), we have 

4 2 Wp 2 ) = ^4 2) (Pi,P2), (3-9) 

Pl2 
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and the second order longitudinal displacement field is then 

*f{k,t) = -\D\t)L {2 \k) 1 (3.10) 

with 

L i2 \k) = J ^^(27r) 3 5g ) (p 12 -fc)gB[ 2) (p 1 ,p 2 )5o(p 1 )5o(p 2 ). (3.11) 
This simple procedure works as long as the displacement field is fully longitudinal in the 

(2) 

Lagrangian frame. However, since p 12 • S T (p 1 ,p 2 ) = due to the transverseness condition, 

(2\ 

eq. (3.5) gives us the longitudinal solution p 12 • S L (p 1 ,p 2 ) only. Thus, even if there were a 
transverse part already at second order, it would not affect the second order density contrast. 
This leads to an important consequence for the density contrast at arbitrary order n: the 
transverse kernel cannot be obtained via the relation (3.7), but it will affect the nth 

order density contrast 5^ if m < n. Thus, (or T m ) has to be constrained at each 

order. 



4 Lagrangian transverse fields 



The Eulerian irrotationality condition states the vanishing of the Eulerian curl of the particle 
motion [9, 23]: 

V x x u(x,t) = 0. (4.1) 

As mentioned in the introduction this implies a restriction to sufficient large scales, where 
we expect that the time evolution should not generate vorticities. Because of the non-trivial 
transformation to Lagrangian coordinates, however, the displacement field must include a 
transverse component as well [13]. Crucially, longitudinal and transverse fields are dynami- 
cally coupled in the Lagrangian picture, 1 even if the transverse fields are zero — the latter is 
the case at lower orders. At higher orders, transverse fields attach the same importance as 
longitudinal fields. 

It is possible to write the irrotationality condition directly in Fourier space, and it leads 
to the general result for n > 2: 



f {n \k) 



d 3 pid 3 £>2 

(27T) 6 



(2vr) 3 5$ (p 12 - k) [p! (p 2 • p 12 ) - p 2 {p 1 ■ p 12 )} /p 



l<i<7, 
i+j=n 



)l 



Z W ( Pl ) + T ( V)l • \L U \p 2 ) + f ij \p 2 ) 



2 

12 



(4.2) 



or, alternatively 



2» 



(*>-£*/ 



^^(2^g>( Pl2 - fc )( Pl x P2 ) 



(2„) 
{ n 



L {j \p 2 ) + t {j \p 2 ] 



1<*<7, 



(4.3) 



1 In SPT the velocity field can be written in longitudinal and transverse fields, and they are decoupled from 
each other. 
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This representation of the original irrotationality condition (4.1) is new, since it is exact at 
any order in perturbation theory and it embeds transverse sources as well; 2 additionally, it 
contains a technical simplification with respect to the commonly used irrotationality condi- 
tion which we shall highlight in the appendix A. With the above expression, the transverse 

~ (n) ~ (n) 

displacement field with time evolution is then (k, t) = —iD n {t) T (k). The perturbation 
vectors in the curly brackets are given by eqs. (2.9) and (2.10), thus the explicit Lagrangian 
formalism is not needed and only the (lower order) results have to be plugged in. 

~ ( 2 ) 

With the use of the Zel'dovich approximation, we immediately obtain for n = 2: T =0 

(2) 

and thus Sj, =0, because of the trivial condition j — i = but also due to symmetry reasons 
in the integrations over p 1 and p 2 - In general, there are vanishing contributions for i = j in 
the summation. Note that in a very general treatment, it is possible to obtain non-vanishing 
transverse solutions already at the second order [25]. We are not considering this case here. 



5 The recursive procedure 

Assuming that eq. (3.7) is valid for arbitrary order n, we obtain the nth order longitudi- 

(s) 

nal displacement field iteratively in terms of Fn and lower order LPT results. Projecting 
out the longitudinal part, i.e., B^f > {p 1 , . . . ,p n ) =Pi... n ■ S™ (Pi> • • • >P n ) > we obtain for the 
longitudinal part of the displacement field 

S f = ^B^\ <=FW-^), (5.1) 

Pl-n 



with 



el" , <*> i ± [^i si ^ + ^ + si ^ 2 ')T . . (5 . 2) 



a=l 



This is our second main result. The last line is the strict consequence of Taylor expanding 
eq. (2.5), and we have introduced the operator On\x} which extracts the symmetric nth 
order part of its argument X, e.g. 

{k ■ S^( Pl ) +k ■ S«( Pl ) k ■ SW(p 2 ,p 3 )} =l{k ■ S«( Pl ) k ■ S®(p 2 ,p 3 ) 

+ two perms. | . (5.3) 

The dependence of the vectors in eq. (5.2) is S^ T = S^ T (p 1 , . . . ,p&), for products of two 

vectors it is p 1 ... k -S { ^ T p 1 ... k -S { [l )T = p 1 ... k -Sf (S)T {p 1 , . . . ,p i )p 1 ... fc -5^ T (p J - +1 , . . . ,pj), with 
k = i + j, and similar for higher order products. 

It is then straightforward to calculate higher order displacement fields. 



6 Example: The third order displacement field 

In this section we demonstrate how the third order displacement field can be obtained and 
expressed in terms of the SPT kernels Fn \ We start with eqs. (5.1) and (5.2) for n = 3, this 

2 Reference [24] gives a Fourier expression for the irrotationality condition as well, however it is non-exact. 



- 6 - 



leads to the longitudinal solution 

4 3) = , B?\ Pl ,p 2 ,p 3 ) = F^\ Pl ,p 2 ,p 3 ) - Ei s \ Pl ,p 2 ,p 3 ) 

Pl23 



with 



+ H + ^P) ( x + ^ir) ^ 2)(p2 ' p3) + tw ° perms j ' (6 - 2) 

where s[ 2) = F 2 W - is given in eq. (3.8). To obtain the longitudinal displacement field 
at third order we use eqs. (2.7) and (2.9). This leads to 

* { ?\k,t) = -iD 3 (t) J d3 P^fP 3 (2vr) 3 ( 5g ) (p 12 3 - fc) Sofa) ftfa) Up 3 ) 

xZ&B® ( Pl ,p 2 ,p 3 ). (6.3) 

^123 

On the other hand, the transverse field (4.2) at third order is 

T ,3 V)=/ ^H^^gW -fc^xfexp,) £<%).£<"(*). (6.4) 

The only thing we have to do is to use the lower order results and and substitute 
the integration limits in the above expression. We then have for the transverse displacement 
field (2.8) at third order 



?( 2 )/„ ^ ^ Pi 'P23 



x S^(p 2) p 3 ) *o(Pi) 5 (P 2 ) SoiPs) , (6-5) 

Pl^23 



and thus we obtain 

* (3 \k,t) = * { * ) (k,t) + * { T ) (k,t). (6.6) 

In general, the use of the recursion relation reduces the work significantly. The final expres- 
sions for higher order displacement fields are surely longer, but the procedure is exactly the 
same compared to the above. 



7 Conclusions 



For the first time, we have formulated an iterative procedure to calculate the Fourier trans- 
form of the Lagrangian displacement field up to arbitrary order in perturbation theory. Our 
procedure is based on the physical assumption that the density contrast agrees in SPT and 
LPT, if the treatment is perturbative and if we restrict our formalism to the initial position 
limit (IPL) (in the IPL, the linear density contrast is evaluated in the vicinity of the initial 
Lagrangian position instead of the evolved Eulerian coordinate; see the thorough discussions 
in [21, 26]). This allows us to relate the LPT series to its counterpart in SPT through 
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the density contrast relation (2.5), and the SPT results are given by the well known SPT 
recursion relation. 

Even for an irrotational motion in the Eulerian frame, the Lagrangian displacement 
field consists not only of longitudinal components but also of transverse components; the 
transverse perturbations affect the nth order density contrast if they are of lower order than 
n. As a consequence, the transverse perturbations cannot be calculated within the density 
contrast relation, but have to be constrained at each order. We have calculated a new 
representation of the irrotationality condition directly in Fourier space, eq. (4.2). This new 
representation has the big advantage that the explicit Lagrangian formalism is not needed 
and only the (lower order) results have to be plugged in. The calculation of the Lagrangian 
transverse fields is straightforward, and so is then the full (i.e., longitudinal and transverse) 
displacement field at any order. 

Some remarks for applications are appropriate here. First of all, the Lagrangian solution 
always contains more non-linear information than the standard one due to the inherent non- 
linearity of the unexpanded relation of the density contrast (2.4). In a future project we shall 
introduce a numerical treatment of the very non-perturbative expression, and we will clarify 
the performance of higher order LPT solutions. Furthermore, the use of our result is not 
restricted to the IPL: One may relax this approximation after the iterative procedure, thus 
effectively readjusting the inherent level of non-linearities (the kernels derived in the IPL are 
still valid). Finally, (higher order) LPT solutions are for example needed for resummation 
techniques of matter poly-spectra (e.g. [27, 28]). 
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A The derivation of equation (4.2) 

The requirement of an irrotational fluid motion yields the Lagrangian constraint equation 
[21]: 

e m^kj ~ e ijk * u ±* ltk = %, n e njk (%j^i,k ~ ^kjj > ( A -!) 

with \I/ = ty(q,r]) being the non-perturbative displacement field. Summation over repeated 
indices is implied, the subscript \f denotes a partial derivative with respect to the Lagrangian 
coordinate qj, and we use the superconformal time dr] = dt/a 2 . Note that SijkL^J = because 

of L^ n > = Vg</>( n ) and due to the symmetry of the second derivatives. Equation (A.l) is the 
strict result of the transformation into Lagrangian coordinates. However, the term in brackets 
on the RHS is essentially redundant, since it is always of higher order. The reason for this 
can be understood, if we rewrite the above equation in a schematic but perturbative way: 

C\ n) = - E , (A.2) 

p+q=n 

where we have defined 

E«w(s*S-*Ss*iS)- (A - 3 » 
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The leading order solution is therefore C\ =0 at any order n. The nth order solution of 
eq. (A.l) is thus 

%^<»> _ £ s m ( £ W + *») - ( £ <«> + T«>) . (A.4) 

' p+q=n J 1 

Since is purely transverse, we can write it in terms of a vector potential: = £kimA^\. 
Then, we Fourier transform the above equation and multiply it with an additional Levi-Civita 
connection. This gives an equation for dT^/dr/. We explicitly evaluate the time derivative 
on the RHS with the use of eq. (2.8), and then integrate (A.4) over superconformal time, 
which yields the final expression (4.2). Note that we have fixed the integration limits in 
eq. (4.2) such that only the fastest growing mode is retained. 
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